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Abstract 

CO . The condition for pressure isotropy is reduced to a recurrence equation 

' with variable, rational coefficients of order three. We prove that this dif- 

, ference equation can be solved in general. Consequently we can find an 

' exact solution to the field equations corresponding to a static spherically 

, symmetric gravitational potential in terms of elementary functions. The 

metric functions, the energy density and the pressure are continuous and 
^H, well behaved which implies that this solution could be used to model the 

interior of a relativistic sphere. The model satisfies a barotropic equation 
of state in general which approximates a polytrope close to the stellar 
centre. 



^ ■ 1 Introduction 



Static solutions of the Einstein field equations for spherically symmetric man- 
ifolds are important in the description of relativistic spheres in astrophysics. 
The models generated may be used to describe highly compact objects where 
the gravitational field is strong as in neutron stars. It is for this reason that 
many investigators use a variety of mathematical techniques to attain exact 
solutions. On of the first models, satisfying all the physical requirements for a 
neutron star, was found by Durgal and Bannerji [1]. Now there exist a number 
of comprehensive collections [2, 3, 4] of static, spherically symmetric solutions 
which provide a useful guide to the literature. It is important to note that 
only a few of these solutions correspond to nonsingular metric functions with a 
physically acceptable energy momentum tensor. 

In this paper we seek a new exact solution to the field equations which can 
be used to describe the interior of a relativistic sphere. We rewrite the Einstein 
equations as a new set of differential equations which facilitates the integration 
process. We choose a cubic form for one of the gravitational potentials, which we 
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believe has not been studied before, which enables us to simplify the condition 
of pressure isotropy. This yields a third order rccTurcncc relation, which we 
manage to solve from first principles. It is then possible to exhibit a new exact 
solution to the Einstein field equations. The curvature and thermodynamical 
variables appear to be well-behaved. We also demonstrate the existence of an 
explicit barotropic equation of state. For small values of the radial coordinate 
close to the stellar core the equation of state approximates a polytrope. We 
believe that a detailed physical analysis of our solution is likely to lead to a 
realistic model for compact objects. 



2 Static spacetimes 

Since our intention is to study relativistic stellar objects it seems reasonable, on 
physical grounds, to assume that spacetime is static and spherically symmetric. 
This is clearly consistent with models utilised to study physical processes in 
compact objects. The generic line element for static, spherically symmetric 
spacetimes is given by 

in Schwarzschild coordinates. For neutral perfect fluids the Einstein fleld equa- 
tions can be written in the form 



-[r{l-e-'')]' = p (2a) 
l_e-2A) + ^e-2A ^ p (2b) 
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^-■^^y< ^y'-^ ^■__y'yl = p (2c) 

for the line element (1) where the energy density p and the pressure p are mea- 
sured relative to the comoving fluid 4-velocity it" = e~^S'Q and primes denote 
differentiation with respect to the radial coordinate r. In the field equations (2) 
we are using units where the coupling constant = 1 and the speed of light 
is c = 1. An equivalent form of the field equations is obtained if we use the 
transformation 

(x) = e^""^'^ ,Z{x) = e-2^(^) ,x = Cr^ (3) 

due to Durgapal and Bannerji [1], where A and C are arbitrary constants. 
Under the transformation (3), the system (2) becomes 

i^-2Z = g (4a, 

4:Zx^y + 2Zx^y+{Zx-Z + l)y = (4c) 

where the overdot denotes differentiation with respect to the variable x. Note 
that (4) is a system of three equations in the four unknowns p,p, y and Z. The 
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advantage of this system lies in the fact that a solution can, upon a suitable 
specification of Z(x), be readily obtained by integrating (4c) which is second 
order and linear in y. 

3 A new series solution 

A large number of exact solutions are known for the system of equations (4) that 
model a relativistic star with no charge. Many of these are listed by Stephani 
et al [2] and Finch and Skea [3]. A comprehensive list of static solutions, 
that satisfy stringent conditions for spherically symmetric perfect fluids, was 
compiled by Delgaty and Lake [4]. The Einstein field equations in the form (4) 
are under-determined. From inspection it is clear that the simplest solutions to 
the system (4) correspond to polynomials forms for Z{x). As far as we are aware 
all exact solutions found previously correspond to forms of the gravitational 
potential Z{x) which are linear or quadratic in the independent variable x. 
Our approach here is to specify the gravitational potential Z{x) and attempt 
to solve (4c) for the potential y. In an attempt to obtain a new solution to the 
system (4) we make the choice 

Z = ax^ + l (5) 

where a is a constant. We suspect that the cubic form (5) has not been consid- 
ered before because the resulting differential equation in the dependent variable 
y is difficult to solve; quadratic forms for Z are listed by Delgaty and Lake [4]. 
The quadratic form for the potential Z is simpler to handle and contains the fa- 
miliar Tolman models. With the specified function Z, the condition of pressure 
isotropy (4c) becomes 

2{ax^ + l)y + Sax'^y + axy = 0. (6) 

The linear second order differential equation (6) is difficult to solve when a ^ 0. 
We have not found a solution for a ^ in standard handbooks of differential 
equations. Software packages such as Mathematica have also not been helpful 
as they generate a solution in terms of hypergeometric functions with complex 
arguments. Wc attempt to find a series solution to (6) using the method of 
Frobcnius. As the point a; = is a regular point of (6), there exist two linearly 
independent solutions of the form of a power series with centre x = 0. We 
therefore can write 

oo 

j/(x) = ^c„x" (7) 

n=0 

where the c„ arc the coefficients of the series. For a legitimate solution we need 
to determine the coefficients c„ explicitly. 
Substituting (7) into (6) yields 

4c2 + (12c3 + aco)x + 4(6c4 + aci)x^ + 

oo 

^{a[2n2 +n+ l]c„ + 2(n + 3)(n + 2)cn+3}ic"^^ = 0. 

n=2 
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For this equation to hold true for all x we require 



4c2 = (8a) 

12c3 + aco = (8b) 

6c4 + aci = (8c) 

a[2n^ + n + l]cn + 2{n + 3){n + 2)cn+3 = 0,n>2. (8d) 

Equation (8d)is a linear recurrence relation with variable, rational coefficients of 
order three. General techniques of solution for difference equations are limited 
to the simplest cases and (8d) does not fall into the known classes. However it 
is possible to solve (8d) from first principles. Equations (8a) and (8d) imply 

C2 = C5 = C8 = ■ ■ ■ = 0. (9) 

From (8b) and (8d) we generate the expressions 

a 1 

2.3- + 3 + 1 1 
~ 22 O 5:2''° 

_a3 2.62 + 6 + 1 2.32 + 3 + l 

~ ~¥ 9:61 8:5:2 

It is clear that the coefficients C3, ce, cg, . . . can all be written in terms of the 
coefficient cq. These coefficients generate a pattern and we can write 

/a\-+ifr 2(3fe)2 + 3fc + l 
C3„^-3 = (-1)"+^ (-) n (3fe + 3)(3fc + 2) ^° ^''^ 

where wc have utilised the conventional symbol Y\ to denote multiplication. 
We can obtain a similar formula for the coefficients C4, C7, cio, . . . From (8c) and 
(8d) we have 

a2.l2 + 1 + 1 

= -2 43 

0^2.42+4 + 12.12 + 1 + 1 

~ 22 TA. 6:3 

_ 0^2.72 + 7 + 12.42 + 4 + 12.12 + 1 + 1 

" ~23 i(x7l ^ i 

The coefficients C4, cy, cio, . . . can all be written in terms of the coefficient c\. 
These coefficients generate a pattern which is clearly of the form 

, ..„+i/a\"+^-rT 2(3fc + l)^ + (3fc + l) + l 

y n (3fc + 4)(3A. + 3) ^^^^ 

where W denotes multiplication. 

Prom (9) we observe that the coefficients 02,05,03,... all vanish. The coef- 
ficients C3, C6, cg, . . . are generated from (10). The coefficients 04,07,010, .. . are 
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generated from (11). Hence the difference equation (8d) has been solved and 
all non-zero coefficients arc expressible in terms of the leading coefficients cq 
and ci- We can write the series (7) as 

y{x) = Co + CiX^ + C3X^ + Cix'^ + C^X^ + CyX^ + CgX^ + CiqX^^ H 



3n+4 



ra=0 / \ n=0 



Co 



n=0 fc=0 



where co and ci are arbitrary constants. Clearly (12) is of the form 

y{x) = coyi(x) + ciy2{x) (13) 



where 



»w ^ (^^E(-o-'(f)"-'n ^'";!r:;g:r' -^i 

are linearly independent solutions of (6). Therefore we have found the general 
solution to the differential equation (6) for the particular gravitational potential 
Z given in (5). The advantage of the solutions in (13) is that they are expressed 
in terms of a series with real arguments unlike the complex arguments given by 
software packages. 



4 Physical models 

Prom (13) and the Einstein field equations (4) we generate the exact solution 

e'^ = } , (14a) 
ax^ + l ^ ' 

^ = -7ax2 (14c) 
^ = 4{ax^ + 1)'^ + ax"^ (14d) 

y 

In the above the quantity y is given by (13), x = Cr'^ and a is a constant. This 
solution has a simple form and is expressed completely in terms of elementary 
functions. The expressions given above have the advantage of simplifying the 
analysis of the physical features of the solution, and will assist in the description 
of a relativistic compact bodies such as neutron stars. 
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Consider a relativistic sphere where < x < b where h = CB? and R is 
the stellar radius. We note that the functions u and A have constant values 
at the centre x = 0. The function p vanishes at the centre. The pressure 
p has a constant value at x = 0. Hence the gravitational potentials and the 
matter variables are finite at the centre. Since y{x) = CQyi{x) + Ciy2{x) is a well 
defined series on the interval [0, h] the quantities v, A, p and p are nonsingular 
and continuous. If a < then the energy density p > 0. The constants cq and 
ci can be chosen such that the pressure p > 0. Consequently the energy density 
and the pressure are positive on the interval [0,6]. At the boundary x = b we 
must have 

g-2A(R) ^ ^^3^6 + 1 ^ 1 _ ^ 

R 

for a sphere of mass M; this ensures that the interior spacetimc matches 
smoothly to the Schwarzschild exterior. For the speed of sound to be less 
than the speed of light we require that 

0<^<1 
dp 

in our units. This inequality will constrain the values of the constants a, cq, 
ci, A and C. Prom this qualitative analysis we believe that the solution found 
can be used as a basis to describe realistic relativistic stars. We believe that 
a detailed physical analysis is likely to lead to realistic models for compact 
objects. 

Our solution has the interesting feature of admitting an explicit barotropic 
equation of state. We observe from (14c) that 



and the variable x can be written in terms of p only. The function y in (13) 
can be expressed in terms of p and the variable x is eliminated. Consequently 
the pressure p in (14d) is expressible in terms of p only, and we can write 

p = p{p). 

Thus the solution in (14) obeys a barotropic equation of state. This highly 
desirable feature is unusual for most exact solutions as pointed out in Stephani 
et al [2]. For small values of x close to the stellar centre we have y ~ cq + cix. 
Then from (14c) and (14d) we have the approximation 

^ (15) 



Therefore for small values of x close to stellar centre (15) implies that we have 
the approximate equation of state 



p (X p 



which is of the form of a polytrope. 
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We point out that the solutions presented in this paper may be extended 
to anisotropic matter. In recent years a number of researchers have proposed 
models corresponding to anisotropic matter where the radial component of the 
pressure differs from the angular component. The physical motivation for the 
analysis of anisotropic matter is that anisotropy affects the critical mass, critical 
surface redshift and stability of highly compact bodies. These investigations are 
contained in the papers [5, 6, 7, 8, 9, 10, 11, 12], among others. It appears that 
anisotropy may be important in fully understanding the gravitational behaviour 
of boson stars and the role of strange matter with densities higher than neutron 
stars. Mak and Harko [11] and Sharma and Mukherjee [13] have observed that 
anisotropy is a crucial ingredient in the analysis of dense stars with strange 
matter. The simple form of our solutions allows for extension to study such 
matter by adapting the energy momentum tensor to include both radial and 
tangential pressures. 
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